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Abstract 

The first-order phase transition in the one-dimensional g-state Potts model with 
long-range interactions decaying with distance as l/r 1+cr has been studied by Monte 
Carlo numerical simulations for < a < 1 and integer values of q > 2. On the basis 
of finite-size scaling analysis of interface free energy AFl, specific heat and Binder's 
fourth order cumulant, we obtain the first-order transition which occurs for a below 
a threshold value o~ c (q). 
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The subject of our study is the one- dimensional (Id) Potts model with ferromag- 
netic long-range (LR) interactions decaying with distance as l/r 1+<J , defined by the 
Hamiltonian 

H = -T T—^6(s hSj ), (1) 
i<j 1 1 J\ 

where J > 0, Sj denotes the g-state Potts variable at site i, S is Kronecker symbol, 
and summation is taken over all pairs in the system. The phase transition at nonzero 
temperature, shown rigorously fij for the Ising (q = 2) case with a < 1 and by 
renormalization group for the continuous n-component models with a < 1 ||, exists 
also in model (1) for o < 1 and q > |3|, [| and goes through a variety of universality 
classes by variation of q and a. Model (1) has been used as a relevant model for 
describing a number of phenomena involving LR interactions, from spin glasses to 
neural networks, but may also be of interest for possible analogies ||, in some cases 
very direct |J, with short-range (SR) models in higher dimensions. 

An important feature of Potts models with SR interactions is the onset of the 
first-order phase transition for q above some threshold value q c (d), which depends on 
dimensionality |7[]. For example, in d = 2 || and d = 4 |J analytical results yield q c 
equal to 4 and 2 respectively, while for d = 3, the approximate methods give a non- 
integer value for q c , slightly lower than 3 HlPl , so that the 3d 3-state Potts model has 



an extremely weak first-order transition, very difficult to detect |Tl| . Such a distinct 
situations created by variation of d and q make this model a canonical example for 
study of various aspects of temperature-driven first-order phase transitions |12| . 

One can expect similar behaviour to occur for the Potts model with LR inter- 
actions, although certain related quantities (e.g. the interface energy) may require 
different interpretation. Generally, this model has been much less explored than its 
SR version, due to non-locality of interactions, which makes difficult the application 
of standard renormalization group techniques in direct space, otherwise appropriate 
for discrete models. A few studies that have been done concern mostly the mean-field 
region and its vicinity by using an e— expansion within Ginzburg-Landau continu- 



ous formalism fl3]|, and the special case of a = 1 B. We have recently proposed 
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the finite-range scaling (FRS) approach |H| suitable for study of the Hamiltonian 
(1), with both a and q arbitrary and continuous 0. This approach is, however, in- 
herently insensitive to discern the first-order transitions and remained inconclusive 
with this respect. The problem was not resolved in other recent works [15|, [16[ on 
the LR Potts model either. 

For this reason we present here the results based on simple numerical simulations 
performed with intention to examine the existence of a first-order transition in this 
model and get a qualitative estimate of its dependence on q and a. 

It has been recently pointed out |I7|, that the temperature-driven first-order 
transitions can be identified from finite-size scaling (FSS) analysis of maxima of the 
energy probability distribution 

Pl(E) = ^L^Um e~ KE , (2) 

where K = J/ksT, Zl(K) is the partition function, Ml{E) is the number of config- 
urations with energy E, and index L denotes the system size. Due to coexistence of 
phases, at the temperature of a first-order phase transition Pl{E) has two maxima, 
corresponding to the two wells in the free energy. The barrier separating them, 
which represents the interface free energy is defined by 

Pl i.E m in) 



AFr. = In 



(3) 



P L (E m , 

where E m i n and E max denote the energies corresponding to the minimum and one of 
the two maxima respectively, the finite-size temperature being adjusted so as to 
make the two maxima equal. For a first-order phase transition AFl should diverge 
with L. In systems with SR interactions it scales like a surface, i.e. ~ while 
in the present case it is expected to scale rather like a volume, i.e. ~ L. 

The calculations were performed on chains of size 100 < L < 400 with periodic 
boundary conditions. We have used the simple Metropolis single-spin-flip algorithm 
with 1 x 10 6 — 3 x 10 6 Monte Carlo (MC) sweeps per spin. The number of necessary 
runs for the precise localization of each of the size-dependent critical temperatures 



Kl has been reduced by applying the Ferrenberg and Swendsen |18| histogram 
method. 
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We have considered integer values of q > 2 in the interval < a < 1. While for 
q = 2 (Ising case), the simulations at T c show only a single maximum in Pl{E) in the 
entire range of a, for higher values of q, the two peaks emerge for o sufficiently low. 
They become more pronounced with increasing q or decreasing a. For illustration, 
in Figure 1 are shown the maxima corresponding to different values of a, taken from 
three typical sets of simulations with fixed q = 3 and L = 400. 

We report here the systematic results for two chosen values, q = 3 and q = 5, in 
the whole interval < o < 1 taken with increment 0.1. 

In Figures 2a and b are summarized the results for the free energy barrier plotted 
as a function of chain size L for q = 3 and q = 5 respectively. The corresponding 
critical temperatures extrapolated to L — > oo, given in table 1, are found to be 
in good agreement with our earlier FRS results [|J, as well as with other known 
approximate results [TB[ |. 



a q Kf c K F e RS q Kf c K F e RS 



0.1 3 0.190 0.136 5 0.262 0.28 

0.2 0.279 0.270 0.333 0.45 

0.3 0.380 0.386 0.492 0.576 

0.4 0.489 0.494 0.637 0.690 

0.5 0.771 0.803 

0.6 0.901 0.920 

0.7 1.019 1.046 



Table 1: Inverse critical temperatures (K^ 1 ), obtained by extrapolation of 
compared to FRS extrapolated values (K FRS ) [Q . 

For both considered values of q, there is a wide range of a, where AFl increases 
with size, indicating the first-order transition. As expected for the LR interactions, 
AFl is proportional to volume rather than surface and depends linearly on L. The 
slope is larger by an order of magnitude for q = 5 in comparison to q = 3, show- 
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ing that the first-order character becomes stronger with increasing q, like for the 
SR interactions. In both cases the slope decreases with increasing a to the point 
where AFl becomes of the order of error bars. Beyond this value, at least for sizes 
considered here, the Pl{E) exhibits a single maximum indicating the onset of the 
second-order phase transition. In present calculation, taken with a rough increment 
of 0.1 in a, this change is observed around a = 0.5 and a = 0.8 for q = 3 and 
q — 5, respectively. These values should be taken with caution and only as a lower 
limit for the threshold value a c between the first- and the second-order transition. 
Namely, the present model allows the continuous approach to the threshold value 
cr c , whereby the first-order transition becomes arbitrarily weak and very difficult to 
detect, comparable to the situation with the 3c? 3-state Potts model with SR inter- 
actions. The above results, however, strongly suggest that a c is considerably larger 
for q = 5 than for q = 3, and that dependence a c (q) should be expected, analogous 
to threshold dependence d c (q) in the SR model. 

Two other energy-related quantities are more conventionally [plj used for de- 
termination of the first-order transition in context of FSS analysis of MC simu- 
lation results: specific heat and Binder's fourth order cumulant H20fl , which both 
can be derived from Pl(E), and expressed in terms of higher energy momenta 
(E n ) L = J2 E E n P L {E) . Specific heat is given by 

Cl = ^( (E 2 ) l - (Ef L ) . (4) 

According to the FSS theory, for second-order transitions its maximum scales as 
(jmax ^ i^a/u^ w h ere a anc i v are fa e critical exponents of the specific heat and the 
correlation length respectively. When the transition is of the first order, it scales as 
a volume, i.e. C™ ax ~ L d . Instead of the Binder fourth cumulant V^f* = 1 — /3, 
we consider here the ratio 



hi) = 5 /;l >J 

< 

r(4) 



For the first-order transitions, lim^oo U L = 1 when T ^ T c , while at T = T c 
lirn^oo = const > 1. For second-order transitions it always tends to one. 
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We present the behaviour of those two quantities on two examples: q = 5, a = 
0.2 and q = 3, o = 0.8, representative of the first- and second-order regimes, 
respectively. 

In Figure 3a and b one can observe two kinds of behaviour of C™ ax in the two 
cases: linear and power law. The fit to the form C™ ax ~ L x in the latter case gives 
the value x = 0.24 for q = 3, a = 0.8. The bare extrapolation error bars for x are 
estimated to be of order of 10%. The hyper-scaling relation with substitution of 
FRS result || v = 1.74 gives ajv = 0.15. The difference can be attributed to the 
general difficulty in extracting the critical exponents from the specific heat, and to 
the fact that the calculated exponent is small and additional correction terms due 
to finite-size gain importance. 

The convergence of the maxima of with size is presented in figure 2 ((c) and 
(d)). The points for the Ising model with o = 0.5 are added as a reference for 
second-order transition behaviour. The case q = 3, a = 0.8 shows, within numerical 
error bars, clear convergence towards 1. The points for q = 5, a = 0.2 converge 
towards a much larger value, which is approximately 2.4 when we take into account 
larger values of L and use the linear extrapolation. 

Thus, the two quantities confirm earlier conclusions based upon the behaviour 
of AF^. However, at present stage, we could not extract from these quantities any 
better precision in determination of cr c (q), so we do not reproduce any systematic 
study for them. 

In summary, by simple numerical calculations, in combination with FSS argu- 
ments, we have shown that the Id LR Potts model for integer q > 2 exhibits the 
first-order phase transition for a below some threshold value a c generally depending 
on q. First-order character becomes weaker with the increase of a, which repre- 
sents a continuous parameter leading from first- to second-order phase transition 
regime. More intensive numerical approach [[2lJ should be needed in future in order 
to determine the threshold value cr c (q). 
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Figure captions: 



Figure 1: 

Maxima of In Pl(E') for q = 3, L = 400 and a = 0.3, 0.4 and 0.5, taken at respec- 
tive values of K L . E' = E/\E (a)\, where the E (a) stay for the zero-temperature 
energies. 

Figure 2: 

MC results of the free energy barrier AFl for sizes L = 100 to 400 for: (a) 
q — 3,<j = 0.1 to 0.4, (b) q — 5, a — 0.1 to 0.7. Notice that the slope is by an 
order of magnitude larger in the case (b). The size of the numerical error bars is 
comparable to or smaller than the size of the points. 

Figure 3: 

MC data for sizes L = 100 to 400: specific heat maxima are plotted for (a) 
q — 3, a — 0.8, (b) q — 5, a — 0.2; ratio maxima are plotted in function 
of inverse size for: (c) q — 3, a — 0.8, (d) q — 5, a — 0.2 (full circles). The 
diamonds correspond to q — 2, a — 0.5 taken as a reference. Lines show the linear 
extrapolations. Notice the common scale on x-axis, but different scales on y-axis. 
The size of the numerical error bars is smaller than the size of the points. 
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Fazni prijelaz prvog reda u Id Pottsovom modelu s dugodoseznim medjudjelovanjem 
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Sazetak 

U jednodimenzionalnom Pottsovom modelu s q stanja i s dugodoseznim medju- 
djelovanjima koja opadaju s udaljenoscu kao l/r 1+cr , Monte Carlo simulacijama je 
promatran fazni prijelaz prvog reda za < a < 1 i cjelobrojne vrijednosti q > 2. Na 
temelju scaling analize slobodne energije medjuplohe, specificne topline i Binderovog 
kumulanta cetvrtog reda, dobivamo prijelaz prvoga reda za a manji od granicne 
vrijednosti cr c (g). 
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Abstract 

The first-order phase transition in the one-dimensional g-state Potts model with 
long-range interactions decaying with distance as l/r 1+<J has been studied by Monte 
Carlo numerical simulations for < a < f and integer values of q > 2. On the basis 
of finite-size scaling analysis of interface free energy AFl, specific heat and Binder's 
fourth order cumulant, we obtain the first-order transition which occurs for a below 
a threshold value cr c (q). 

Physics Abstracts classification number: 05.50.+q, 64.60.Cn 



The subject of our study is the one- dimensional (Id) Potts model with ferromag- 
netic long-range (LR) interactions decaying with distance as l/r 1+<J , defined by the 
Hamiltonian 



where J > 0, denotes the g-state Potts variable at site z, 6 is Kronecker symbol, 
and summation is taken over all pairs in the system. The phase transition at nonzero 
temperature, shown rigorously [1] for the Ising (q = 2) case with a < 1 and by 
renormalization group for the continuous n-component models with a < 1 [2], exists 
also in model (1) for a < 1 and q > [3, 4] and goes through a variety of universality 
classes by variation of q and a. Model (1) has been used as a relevant model for 
describing a number of phenomena involving LR interactions, from spin glasses to 
neural networks, but may also be of interest for possible analogies [5], in some cases 
very direct [6], with short-range (SR) models in higher dimensions. 

An important feature of Potts models with SR interactions is the onset of the 
first-order phase transition for q above some threshold value q c (d) } which depends on 
dimensionality [7]. For example, in d = 2 [8] and d = 4 [9] analytical results yield q c 
equal to 4 and 2 respectively, while for d = 3, the approximate methods give a non- 
integer value for q C} slightly lower than 3 [10], so that the 3 d 3-state Potts model has 
an extremely weak first-order transition, very difficult to detect [11]. Such a distinct 
situations created by variation of d and q make this model a canonical example for 
study of various aspects of temperature-driven first-order phase transitions [12]. 

One can expect similar behaviour to occur for the Potts model with LR inter- 
actions, although certain related quantities (e.g. the interface energy) may require 
different interpretation. Generally, this model has been much less explored than its 
SR version, due to non-locality of interactions, which makes difficult the application 
of standard renormalization group techniques in direct space, otherwise appropriate 
for discrete models. A few studies that have been done concern mostly the mean-field 
region and its vicinity by using an e— expansion within Ginzburg-Landau continu- 
ous formalism [13], and the special case of a = 1 [3]. We have recently proposed 




J 
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the finite-range scaling (FRS) approach [14] suitable for study of the Hamiltonian 
(1), with both a and q arbitrary and continuous [4]. This approach is, however, in- 
herently insensitive to discern the first-order transitions and remained inconclusive 
with this respect. The problem was not resolved in other recent works [15, 16] on 
the LR Potts model either. 

For this reason we present here the results based on simple numerical simulations 
performed with intention to examine the existence of a first-order transition in this 
model and get a qualitative estimate of its dependence on q and a. 

It has been recently pointed out [17], that the temperature-driven first-order 
transitions can be identified from finite-size scaling (FSS) analysis of maxima of the 
energy probability distribution 

Pl(E) = ^-fdE) e~ KE , (2) 

where K = J/IcbT, Zl(K) is the partition function, Nl(E) is the number of config- 
urations with energy E, and index L denotes the system size. Due to coexistence of 
phases, at the temperature of a first-order phase transition Pl(E) has two maxima, 
corresponding to the two wells in the free energy. The barrier separating them, 
which represents the interface free energy is defined by 



AFr, = In 



K L 



(3) 



where £ m ,„ and E max denote the energies corresponding to the minimum and one of 
the two maxima respectively, the finite-size temperature Kl being adjusted so as to 
make the two maxima equal. For a first-order phase transition AFl should diverge 
with L. In systems with SR interactions it scales like a surface, i.e. ~ L d ~\ while 
in the present case it is expected to scale rather like a volume, i.e. ~ L. 

The calculations were performed on chains of size 100 < L < 400 with periodic 
boundary conditions. We have used the simple Metropolis single-spin-flip algorithm 
with 1 X 10 6 — 3 X 10 6 Monte Carlo (MC) sweeps per spin. The number of necessary 
runs for the precise localization of each of the size-dependent critical temperatures 
Kl has been reduced by applying the Ferrenberg and Swendsen [18] histogram 
method. 
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We have considered integer values of q > 2 in the interval < a < 1. While for 
q = 2 (Ising case), the simulations at T c show only a single maximum in Pl(E) in the 
entire range of <r, for higher values of q, the two peaks emerge for a sufficiently low. 
They become more pronounced with increasing q or decreasing a. For illustration, 
in Figure 1 are shown the maxima corresponding to different values of <r, taken from 
three typical sets of simulations with fixed q = 3 and L = 400. 

We report here the systematic results for two chosen values, q = 3 and q = 5, in 
the whole interval < a < f taken with increment O.f . 

In Figures 2a and b are summarized the results for the free energy barrier plotted 
as a function of chain size L for q = 3 and q = 5 respectively. The corresponding 
critical temperatures extrapolated to L — > oo, given in table 1, are found to be 
in good agreement with our earlier FRS results [4], as well as with other known 
approximate results [16]. 



a q Kf c K F e RS q Kf c K FRS 



0.1 3 0.190 0.136 5 0.262 0.28 

0.2 0.279 0.270 0.333 0.45 

0.3 0.380 0.386 0.492 0.576 

0.4 0.489 0.494 0.637 0.690 

0.5 0.771 0.803 

0.6 0.901 0.920 

0.7 1.019 1.046 



Table 1: Inverse critical temperatures (K^ c ), obtained by extrapolation of Kl 
compared to FRS extrapolated values (K FRS ) [4]. 

For both considered values of q } there is a wide range of <r, where AFl increases 
with size, indicating the first-order transition. As expected for the LR interactions, 
AFl is proportional to volume rather than surface and depends linearly on L. The 
slope is larger by an order of magnitude for q = 5 in comparison to q = 3, show- 
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ing that the first-order character becomes stronger with increasing q, like for the 
SR interactions. In both cases the slope decreases with increasing a to the point 
where AFl becomes of the order of error bars. Beyond this value, at least for sizes 
considered here, the Pl(E) exhibits a single maximum indicating the onset of the 
second-order phase transition. In present calculation, taken with a rough increment 
of 0.1 in cr, this change is observed around a = 0.5 and a = 0.8 for q = 3 and 
<? = 5, respectively. These values should be taken with caution and only as a lower 
limit for the threshold value o c between the first- and the second-order transition. 
Namely, the present model allows the continuous approach to the threshold value 
cr c , whereby the first-order transition becomes arbitrarily weak and very difficult to 
detect, comparable to the situation with the 3 d 3-state Potts model with SR inter- 
actions. The above results, however, strongly suggest that o c is considerably larger 
for q = 5 than for q = 3, and that dependence cr c (q) should be expected, analogous 
to threshold dependence d c (q) in the SR model. 

Two other energy-related quantities are more conventionally [19] used for de- 
termination of the first-order transition in context of FSS analysis of MC simu- 
lation results: specific heat and Binder's fourth order cumulant [20], which both 
can be derived from Pl(E) } and expressed in terms of higher energy momenta 
(E n ) L = Y.E Enp L{E) . Specific heat is given by 

Cl = ( (E 2 ) L - (E)l ) . (4) 

According to the FSS theory, for second-order transitions its maximum scales as 

v are the critical exponents of the specific heat and the 
correlation length respectively. When the transition is of the first order, it scales as 
a volume, i.e. C^ ax ~ L d . Instead of the Binder fourth cumulant = 1 — t/j^/3, 
we consider here the ratio 



(4) _ < E 4 > L , , 



For the first-order transitions, limj^oo = 1 when T ^ T c , while at T = T c 
limj^oo Uj 4 ^ = const > 1. For second-order transitions it always tends to one. 



5 



We present the behaviour of those two quantities on two examples: q = 5, a = 
0.2 and q = 3, a = 0.8, representative of the first- and second-order regimes, 
respectively. 

In Figure 3a and b one can observe two kinds of behaviour of C™ ax in the two 
cases: linear and power law. The fit to the form C™ ax ~ L x in the latter case gives 
the value x = 0.24 for q = 3, a = 0.8. The bare extrapolation error bars for x are 
estimated to be of order of 10%. The hyper-scaling relation with substitution of 
FRS result [4] v = 1.74 gives ajv = 0.15. The difference can be attributed to the 
general difficulty in extracting the critical exponents from the specific heat, and to 
the fact that the calculated exponent is small and additional correction terms due 
to finite-size gain importance. 

The convergence of the maxima of with size is presented in figure 2((c) and 
(d)). The points for the Ising model with a = 0.5 are added as a reference for 
second-order transition behaviour. The case q = 3, a = 0.8 shows, within numerical 
error bars, clear convergence towards 1. The points for q = 5,<r = 0.2 converge 
towards a much larger value, which is approximately 2.4 when we take into account 
larger values of L and use the linear extrapolation. 

Thus, the two quantities confirm earlier conclusions based upon the behaviour 
of AFl. However, at present stage, we could not extract from these quantities any 
better precision in determination of cr c (q) } so we do not reproduce any systematic 
study for them. 

In summary, by simple numerical calculations, in combination with FSS argu- 
ments, we have shown that the Id LR Potts model for integer q > 2 exhibits the 
first-order phase transition for a below some threshold value o c generally depending 
on q. First-order character becomes weaker with the increase of <r, which repre- 
sents a continuous parameter leading from first- to second-order phase transition 
regime. More intensive numerical approach [21] should be needed in future in order 
to determine the threshold value cr c (q). 
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Figure captions: 



Figure 1: 

Maxima of In Pl(E') for q = 3, L = 400 and a = 0.3, 0.4 and 0.5, taken at respec- 
tive values of Kl. E' = E /\E (cr)\ } where the E (cr) stay for the zero-temperature 
energies. 

Figure 2: 

MC results of the free energy barrier AFl for sizes L = 100 to 400 for: (a) 
q = 3,<T = 0.1 to 0.4, (b) q = 5,<r = 0.1 to 0.7. Notice that the slope is by an 
order of magnitude larger in the case (b). The size of the numerical error bars is 
comparable to or smaller than the size of the points. 

Figure 3: 

MC data for sizes L = 100 to 400: specific heat maxima are plotted for (a) 
q = 3, a = 0.8, (b) q = 5, cr = 0.2; ratio maxima are plotted in function 

of inverse size for: (c) q = 3,<r = 0.8, (d) q = 5,<r = 0.2 (full circles). The 
diamonds correspond to q = 2, a = 0.5 taken as a reference. Lines show the linear 
extrapolations. Notice the common scale on x-axis, but different scales on y-axis. 
The size of the numerical error bars is smaller than the size of the points. 
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Sazetak 

U jednodimenzionalnom Pottsovom modelu s q stanja i s dugodoseznim medudje- 
lovanjima koja opadaju s udaljenoscu kao l/r 1+<J , Monte Carlo simulacijama je 
promatran fazni prijelaz prvog reda za < a < 1 i cjelobrojne vrijednosti q > 2. Na 
temelju scaling analize slobodne energije meduplohe, specificne topline i Binderovog 
kumulanta cetvrtog reda, dobivamo prijelaz prvoga reda za a manji od granicne 
vrijednosti cr c (q). 



